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Here, we present theoretical analysis of electron transport in polyaniline based (PANi) nanofibers 
assuming the metalic state of the material. To build up this theory we treat conducting polymers 
as a special kind of granular metals, and we apply the quantum theory of conduction in mesoscopic 
systems to describe the transport between metallic-like granules. Our results show that the con- 
cept of resonance electron tunneling as the predominating mechanism providing charge transport 
between the grains is supported with recent experiments on the electrical characterization of single 
PANi nanofibers. By contacting the proposed theory with the experimental data we estimate some 
important parameters characterizing the electron transport in these materials. Also, we discuss 
the origin of rectifying features observed in current- volt age characteristics of fibers with varying 
cross-sectional areas. 

PACS numbers: 72.80.-r 



Starting from their discovery, conducting polymers un- 
dergo intense research. These materials are significant 
mostly due to various possible applications in fabrica- 
tion of nanodevices, for polymer based devices should 
have advantages of low cost and flexible, controlled chem- 
istry 0. Also, there are some unsolved problems con- 
cerning physical nature of charge transfer mechanisms 
in conducting polymers making them interesting sub- 
jects to fundamental research. Recently, we reported 
of fabrication of doped polyaniline/polyethylene oxides 
(PANi/PEO) nanofibers using electrospinning technique 
[2J, LI ■ The details of the electrospinning setup used for 
the fibers fabrication are described elsewhere 0, • 

Single fiber electrical characterization was carried out 
at room temperatures (T ~ 300 K). It was revealed in 
the experiments that current- voltage curves for conduct- 
ing PANi fibers could be non-ohmic. Field-effect transis- 
tor (FET) saturation behavior was observed in the device 
made out of two nanofibers with the diameters 300 nm 
and 120 nm, respectively, placed across two gold elec- 
trodes j2j on a doped Si/SiC>2 wafer in a FET config- 
uration. As well, nonlinear current-voltage curves were 
obtained in trasport experiments on thinner nanofibers 
. Some characteristics appeared to be asymmetric with 
respect to the reversal of the applied bias voltage polar- 
ity, especially those for fibers of varying diameter. We 
do believe that these experimental results contain impor- 
tant information concerning physical mechanisms of the 
charge transport in conducting polymers. In the present 
work we concentrate on the theoretical analysis of the 
previous experimental results to recover this information. 

Conducting chemically doped polymers are known to 
be very inhomogeneous. In some regions polymer chains 
are disordered, forming amorphous poorly conducting 
substance. In other regions the polymer chains are or- 
dered and densely packed @, @ . These regions behave as 
metallic-like grains embedded in a disordered amorphous 



environment. This gives grounds to apply the model of a 
granular metal to describe such materials. The frac- 
tion of metallic-like islands in bulk polymers varies de- 
pending on the details of synthesis process. Presumably, 
in actual samples the metallic grains always remain sepa- 
rated by poorly conducting disordered regions and do not 
have direct contacts. Within this model electronic states 
are assumed to be delocalized over the grains, therefore 
electrons behave as conduction electrons in conventional 
metals. Their motion inside the grains is diffusive with 
the diffusion coefficient D — \v 2 f t [vf is the Fermi ve- 
locity, and r is the scattering time). 

Recently, Prigodin and Epstein (PE) suggested that 
the grain-to-grain transport mostly occurs due to the res- 
onance tunneling of electrons. The latter is provided with 
intermediate states on the chains connecting the grains 
||. In the present work we further develop the PE ap- 
proach to the charge transport in conducting polymers, 
and we show that it brings results consistent with those 
obtained in transport experiments on individual PANi 
nanofibers 0>I3- 

Within this approach we immediately see a similarity 
in electron transport mechanisms in conductive polymers 
and those in molecular wires connecting metal leads 
In both cases the transport is mostly provided by inter- 
vening " bridges" giving rise to intermediate states for the 
electron tunneling. Employing the granular metal model 
we can treat metallic islands as " leads" connnected with 
single-site bridges. This simple structure of the bridges is 
proposed in PE theory. We also explain the issue below. 
A systematic theory of conductance in quantum molec- 
ular wires is developed during the last two decades. In 
further analysis we employ this theory. Accordingly, we 
write an electric tunneling current flowing between two 
crystalline domains Ii g in the form |l0| : 
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I ig = -r dET eff (E) ME) - f 2 (E) . (1) 
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Here, fi, 2 {E) are Fermi functions taken with different 
contact chemical potentials for the grains: ji\ and /i 2 , 
respectively. The chemical potentials differ due to the 
voltage SW applied between the grains [TlT |: 

m =E F + (1- r))e5W\ fi 2 = E F - rjeSW. (2) 

The dimensionless parameter r\ characterizes how the 
voltage is divided between the grains, Ef is the equi- 
librium Fermi energy of the system, and T e ff(E) is the 
electron transmission function. 

In general, the electron transmission function T e f f (E) 
includes a term T(E) describing coherent transport by 
electron tunneling along with inelastic contribution orig- 
inating from motions in the medium between the grains. 
These motions, especially nuclear motions in the reso- 
nance chain, could cause the electronic phase breaking 
effect. In our case "the bridge" is a single electron level 
situated between the grains. 

For the considered problem the electron transport is a 
combination of tunneling through two barriers (the first 
one separates the left metallic domain from the interme- 
diate state in the middle of the resonance chain, and the 
second separates this state from the right grain, suppos- 
ing the transport from the left to the right) affected with 
the inelastic scattering at the bridge, as shown below. 
The barriers are represented by the squares, and the tri- 
angle in between imitates a scatterer coupling the bridge 
to a dissipative electron reservoir. 
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An electron could be injected into this system, and/or 
leave from there via four channels indicated in this 
schematic drawing. Incoming particle fluxes (Ji) are re- 
lated to those outgoing from the system (Jj) by means 
of the transmission matrix T 0, ^3 : 



Jj — TjiJi. 



1 < i, 3 < 4. 



(3) 



Off diagonal matrix elements Tji (E) are probabilities for 
the electron to be transmitted from the channel i to the 
channel j, whereas diagonal matrix elements Ta(E) are 
probabilities for its reflection back to the channel i 
To provide the charge conservation, the net particle flux 
in the channels connecting the system with the reservoir, 
must be zero. So, we have: 



h + Ji - J3 - 4 = 0. 



(4) 



We have no grounds to believe that there is difference 
between the channels 3 and 4, therefore we assume J3 = 



J 4 . Also, in the following calculations we put J2 = 0, 
presuming no income flux in the channel 2. 

Now, we employ (3) to express outgoing fluxes in terms 
of incoming ones, and we substitute the results into the 
equation (4). This gives: 



J3 — Ja 



Ji 



K\(E) 
2-R{E)' 



(5) 



Here. 



K 1 (E)=T 31 +T 41 ; 



R{E) = T 33 + T 44 + T 43 + T- 
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The transmission function T e ff(E) relates the particle 
flux outgoing from the channel 2 to that incoming to the 
channel 1, namely: 



J 2 - T eff(E)Ji- 
So, using Eqs.(3),(5) we obtain: 

K x {E)Kz{E) 



Teff = T{E) 



1 - R{E) 



where 



T{E)=T 21 ; K 2 (E)+T 23 +Tr 
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(7) 



(8) 



(9) 



The functions T(E), Ki$(E) and R{E) are expressed in 
terms of the matrix elements of the scattering matrix S 
relating outgoing wave amplitudes b[, b' 2 ,a 3 , a 4 to the in- 
cident ones 61, 62, 03, »4 : Tij — \Sq\ 2 - In the considered 
case the S matrix takes the form [Til IHl : 



a 2 r 2 



s = z- 



at\t 2 f3t\ a(jt\r 2 \ 

atit 2 r 2 + a 2 r 1 a(3rit 2 (3t 2 
(3t\ a(3r\t 2 P 2 ri ar\r 2 — a 
\ a(3tir 2 [3t 2 arir 2 — a f3 2 r 2 J 

(10) 

where Z = 1 — o?r\r 2 \ a = \/l — e; (5 = y/e; r\ t2 and t\ t2 
are the amplitude transmission and reflection coefficients 
for the barriers ( | i 1 ,2 1 2 + |n.,2| 2 = 1); and the parameter 
e characterizes the dephasing strength. This parameter 
takes values within the range [0,1], e — correspond- 
ing to the completely coherent, and e = 1 to the fully 
incoherent transport, respectively. 

When the bridge is detached from the dephasing reser- 
voir T e ff(E) — T(E). In this case we can employ a simple 
analytical expression for the electron transmission func- 
tion nnnniii: 



T(E) = 4Ai(£;)A 2 (i;)|G(£;)| s 



where: 

Ai, 2 (£) 



V 2 

v l,2 



E - ei 2 - o-i, 



(11) 



(12) 
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Self-energy terms £1,2 appear due to the coupling of the 
metallic grains ("leads") to the bridge, and Vj.,2 are the 
corresponding coupling strengths. The form of these 
terms originates from a frequently used model where 
the grains are simulated with semiinfinite homogeneous 
chains. The parameters <j\.i represent self-energy cor- 
rections of the chains. In further analysis we estimate 
Si. 2 at the energies close to the Fermi energy Ep, which 
brings the result: £^2 « —iV? 2/71,2 Here, 71,2 is 

the coupling strength in the chains imitating the metallic 
grains. 

The retarded Green's function for a single site bridge 
could be approximated as follows 



G{E) 



E-E-i 



(13) 



where E\ is the site energy. The width of the resonance 
level between the grains is described with the parameter 

r = Ai + a 2 . 

Comparing the expressions (8) and (11) in the limit e = 
we arrive at the following expressions for the tunneling 
parameters 5\^{E) : 



6 1<2 (E)=tl 2 (E) 



2Ai 



V(£-£i) 2 +T2 



(14) 



Using this result we easily derive the general expression 
for the electron transmission function: 



Teff(E) 



g(£)(l + a 2 )[ g (£)(l + a 2 ) 
[g(E)(l- « 2 ) + l + c 



(15) 



where g{E) = ti{E)t 2 {E). 

In further calculations we keep in mind that the trans- 
mission coefficient for the resonance tunneling is deter- 
mined with the probability of finding the resonance state. 
The latter is estimated as T ~ exp(— L/£) where L is the 
average distance between the interacting grains, and £ 
is the localization length for electrons |8j. It takes val- 
ues small compared to unity but much greater than the 
transmission probability for sequental hoppings along the 
chains Th ~ exp(— 2L/£) [i(J. The probability for exis- 
tence of a resonance state at a certain chain is rather low, 
so only a few out of the whole set of the chains connecting 
two grains are participating in the process of intergrain 
electron transport. Therefore one could assume that any 
two metallic domains are connected with a single chain 
providing an intermediate state for the resonance tun- 
neling. All remaining chains could be neglected for they 
poorly contribute to the transport compared to the res- 
onance chain. 

Realistic PANi nanofibers prepared via electrospinning 
have diameters in the range ~ 20 -7- lOOnm and lengths of 
the order of several microns. This is significantly greater 
than typical sizes of metal-like grains and intergrain dis- 
tances, which take on values ~ 5nm |2j. So, we treat a 
nanofiber as a set of parallel channels for the tunneling 



current, any channel being a sequence of metallic do- 
mains connected with resonance chains. The total cur- 
rent in the fiber is the sum of the contributions from 
all these channels. We assume these contributions to be 
equal. Also, the voltage W applied across the whole fiber 
is a sum of contributions SW from sequental pairs of 
grains along a single channel. We assume as the fist ap- 
proximation SW ~ WL/Lf where Lf is the fiber length. 
Since we presume the current to be conserved in any sin- 
gle channel, we can identify the intergrain current Ii g 
given by (1) with the current in the channel. Using the 
above approximation for SW we express the current in 
the channel in terms of the net voltage W instead of SW. 

For certain doping and crystallinity rates the number 
of channels is proportional to the fiber cross-sectional 
area. In reality, dedoping of the fibers due to the con- 
tact of their surfaces with atmospheric gases affects the 
number of the working channels n. This occurs because 
the electron localization length £ decreases in the de- 
doped surface layers, bringing a crucial fall in the conduc- 
tion through the channels which run there. As a result 
these channels become insulating. The relative number 
of the insulating channels is greater in thinner nanofibers, 
therefore their conduction (the latter is proportional to 
the number of working channels) is much stronger re- 
duced due to the dedoping processes, than the conduc- 
tion of thicker samples. 

Conducting samples studied in the experiments of Q 
included 70nm diameter nanofiber (Sample 1), and a pair 
of 18nm and 2hnm diameter fibers connected in paral- 
lel (Sample 2). PANi nanofibers whose diameters were 
smaller than 15nm appeared to be insulating due to de- 
doping, as reported in [3j. Therefore we assume that 
outer layers whose thickness is ~ 8nm are insulating in 
all fibers used in the experiments. Accepting the value 
5 nm to estimate both average grain size and intergrain 
distance and keeping in mind the effect of dedoping, we 
find that the 70 nm fiber in the experiments of |3J] could 
include about 30 -j- 40 conducting channels and we could 
not expect more than 2 working channels in the pair of 
fibers S2. This gives the ratio of conductions ~ 15 — 20 
which is greater than the ratio of cross-sectional areas of 
the samples {A\/A 2 ~ 5). The difference originates from 
the stronger dedoping of thinner fibers. 

To calculate the current in the whole nanofiber I we 
multiply the channel in a single current by the number 
of working channels / = nli g . In calculations we con- 
sider I Vi ; 2 1 << |ti,2 I because the coupling strength of 
the grains to the intermediate site in between is propor- 
tional to a small factor exp(— L/£). As for the magnitude 
of the coupling strengths 71,2, we can roughly estimate 
them as leV. This value is typical (in order) for covalent 
bonds. To provide a good match of our theory with the 
experimental data we employ the least square procedure. 
This brings Vi,2 = S.OmeV for the IQnm fiber (Sample 
1) for n = 35, and 14,2 = 2.4meU (n = 2) for the pair of 
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FIG. 1: Calculated (dash-dot lines) and experimental 
(solid lines) current (nA) - voltage (V) characteristics for 
PANi/PEO nanohbers (T = 300 K). Experimental charac- 
teristics were reported in y]. 



thinner hbers (Sample 2). 

The value of the dephasing parameter e is determined 
with the interplay between the characteristics of electron- 
phonon interaction and those of the electron coupling to 
the leads. When the latter is weak, as in the consid- 
ered case, e should depend on energy E exhibiting peaks 
which correspond to vibrational modes in the resonance 
chain. This was shown in recent papers 0, 0] bas- 
ing on the nonequilibrium Green's function formalism. 
These peaks should be revealed in the transmission func- 
tion T e ff(E) bringing step-like features in the current- 
voltage characteristics. The absence of such features in 
the experimental curves reported in the work ,3j gives 
grounds to believe that in these experiments the rele- 
vant energy range does not include phonon harmonics 
associated with the vibrational modes. Therefore, in our 
calculations we treat e as a constant. At a distance from 
the peaks, the value of e is mostly determined with the 
strength of stochastic motions in the resonance chain, so 
this value depends on temperature. Basing on the results 
of some earlier works (see e.g. |l3|) concerning electron 
transport through molecular bridges at room tempera- 
ture, we assume e = 0.05 in our calculations. 

Our results are shown in the Fig.l. The curves SI, 

52 present current-voltage characteristics for thick 70nm 
nanofiber (SI), and the pair of thin fibers connected in 
parallel (S2). In plotting the strongly asymmetric curve 

53 we put n = 1, and we obtained Vi,2 = 3.2meV. The 
division parameter value is taken 0.35 for the curve SI. 
This value is chosen to reproduce slight asymmetry of 
the corresponding experimental curve. On the contrary, 
to plot the curve S3 we assume r\ = 0.02 which repre- 
sents very asymmetric voltage division along the reso- 
nance chain. 

The asymmetric curve S3 looks like experimental volt- 
ampere characteristic reported in Q for the fiber whose 
diameter was sharply changing from 70 nm to 20 nm near 
the middle. The latter did show asymmetric strongly 



rectifying volt-ampere curves. We conjecture that sud- 
den change in the fiber diameter produces significant 
changes in the voltage distribution along that segment of 
the fiber where its diameter changes, affecting the volt- 
age division between the adjacent metallic grains located 
there, so that the major part of the voltage fits the thin- 
ner segment. Assuming for simplicity that the change 
in the fiber diameter happens along the segment whose 
length has the same order as intergrain distances, we can 
roughly estimate the voltage division comparing cross- 
sectional areas of the segments. For the fibers used in 
the experiments of j^] this ratio is 12.5. Recalling the 
effects of dedoping on both segments we estimate that 
voltage division between the segments could be ~ 25 : 1 
which gives for the division parameter at the interface 
separating the segments estimation: r\ ~ 0.04. This is 
reasonably close to the value r\ = 0.02 used in plotting 
the curve S3. 

Presented theoretical volt-ampere characteristics rea- 
sonably match the curves obtained in the experiments. 
This justifies the accepted concept of the electron trans- 
port in the polymeric granular metals. Also, basing 
on our results we can estimate values of the coupling 
strengths V\^- We conclude that the coupling strength 
of metallic domains to the intermediate state in PANi 
fibers could be estimated as 2.0 -j- 3.5 me V. 

Finally, in this paper we apply theory of electron trans- 
port through molecular wires to analyze transport char- 
acteristics obtained in experiments on PANi nanofibers. 
Demonstrated agreement between the presented theoret- 
ical results and experimental evidence proves that elec- 
tron tunneling between metallic-like grains through inter- 
mediate states at the resonance chains really plays a sig- 
nificant part in the transport in the conducting polymers 
provided that they are in the metallic state. This is a 
novel result which helps to achieve better understanding 
of physics of metallic state of conducting polymers. Fur- 
ther and more thorough analysis based on the concepts 
used here promises to bring more information concern- 
ing physical mechanisms and characterictics of transport 
processes in these substances. 
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